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Abstract 


The 1 eduction of helicopter vibration remains a major challenge to the rotorcraft designer 
The present study addresses the problem of foimulation and design of a multi-input- 
multi output closed loop vibiation control scheme, of a 3 D flexible fuselage model based 
on the concept of ACSR(Active Control Of Structural Response) The sensoi locations 
for vibration measurement are chosen from a optimal set of locations measuring high 
levels of vibration Initially the dynamic model of the full order system is decomposed 
into flexible and rigid body mode subsystems A reduced order model for the flexible 
mode system is formulated using balanced realisation based-order reduction technique 
Conti ollei IS then designed for this reduced model using factorisation theory and internal 
model principle for disturbance rejection The simulation study showed good vibration 
reduction in the fuselage and gearbox Increasing the number of sensors seemed to reduce 
the magnitude of control forces required for vibration minimisation The choice of low 
vibratory level node location for sensor placement gave unsatisfactoiy results 

It IS obseived that sensor locations have a significant influence on the level of vibration 


reduction m fuselage structure 
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Chapter 1 


INTRODUCTION 

1 1 Introduction To The Vibration Control Problem 

In helicopters, controlling vibrations to acceptable limits is one of the key problems facing 
the rotorcraft designer The time dependant aerodynamic and inertia loads acting on the 
rotor blades during forward flight cause vibration in helicopters The harsh vibrational 
environment results in discomfort to pilot and crew and fatigue damage to the structure 
In addition, vibration contributes to increased maintenance cost, difficulty in reading in 
struments and reduced effectiveness of weapons system A comprehensive review on the 
sources of vibration and description of the sources of vibration and description of the 
methods for reducing vibration levels are given m [14] and [18] 

The increasing demands for passenger transportation, high speed and high perfor 
mance, improved manoevrabihty coupled with the need to improve system reliability 
and reduce maintenance costs has resulted m stringent requirements for vibrations For 
present day helicopters, the general requirement is to have a maximum vibration level of 
0 Ig in the fuselage With the adoption of stringent requirements (as given in [4]), it will 
become necessary m future to reduce vibratory levels below 0 05g or even lower 

The strategies employed for vibration control in helicopters can be broadly classified 
as passive or active control The passive approach utilizes vibration absorbers and vibra 
tion isolation devices Generally, their incorporation leads to increase in rotor weight and 
aerodynamic drag Another passive approach is, careful structural dynamic design using 
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structural optimisation, aimed at minimising vibration in forward flight The inability of 
the passive methods and devices to adapt to the wide operating conditions, represents a 
severe limitation 

The active control methodologies include HHC (Higher harmonic control) IBC (Individual 
blade control), ACF(Actively controlled flap) and ACSR(Active control of structural 
response) which has been described in detail in [9] The higher harmonic pitch control 
of the blade provides a substantial reduction of vibratory levels in the fuselage but cause 
rotor blade stall at high speeds and require high power and low airworthiness criteria 
ACF has been able to eliminate some of these drawbacks 

The ACSR scheme is based on the idea that in a linear system one can superimpose 
two independent response quantities such that total response is zero See e g [3] and [19] 

Vibration reduction is acheived by exciting the fuselage with actuators located at the 
gearbox support mountings such that the sum of the response of the airframe due to 
rotor loads and external excitation is reduced to a minimum Refer [3], [20] Ground and 
flight tests on the Westland 30 helicopter with the ACSR system, provided remarkable 
vibration reduction throughout the whole flight envelope explored See [25] The major 
advantages of ACSR are 

• Ability to minimise vibration at specific fuselage location 

• Simplicity and minimal impact on air worthiness requirements because vibration 
control IS implemented entirely in the non rotating system 

• Low power requirements 

• It can operate at a multiplicity of frequencies 

The analytical study on prediction and reduction of vibration in helicopter requires 
a realistic fuselage model A 3 D finite element model of a helicopter fuselage, shown in 
fig (3 2) developed by Mangalick, et al, [23] has been used in the present study The struc 
tural dynamic characteristics of this model resembles those of a realistic helicopter The 
coupled gearbox fuselage dynamic model is given in fig (4 1) 
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The helicopter fuselage is described by 64 nodes (with 384 degrees of freedom) as 
shown in fig (3 2) Assuming that the rotor system consists of four blades, the vibratory 
hub loads will have a nondimensional excitation frequency of 4/rev For the fuselage 
model the nondimensional natural frequency of the 20th flexible mode is 6 4064 Ref[2l| 
which is 50% more than the excitation frequency (4/rev) of the hub loads Therefore 
the first 20 modes of the helicopter fuselage are considered in the vibration analysis 
Since the vibratory level in the vertical direction(z) in the fuselage is more predominant 
without loss of generality it is assumed that the sensors measure only the z component 
of the fuselage vibration The vibration control is performed only to minimise these z 
component of vibration 

The total number of degrees of freedom is 26 with 3 rigid body modes (pitch heave, roll) 
of gear box, 3 rigid body modes (pitch, heave, roll) of fuselage and 20 flexible modes uf 
fuselage In state space it is described by a 52 order system For the present study the 
vibratory response of the fuselage for only 4/rev disturbance excitation is considered 


1 2 Objectives Of The Present Study 

The mam problems addressed in this thesis are 

• Formulation of a reduced order model, satisfactorily approximating the frequency 
response of the original system 

• Formulation of a MIMO closed loop control scheme for vibration minimization, at 
the prescribed disturbance frequency (4/rev), using reduced order model 

• Analyse the effectiveness of the closed loop control scheme in reducing the vibration 
in the fuselage, by incorporating in the full order model 

The order reduction and closed loop control scheme are acheived by formulating a input 
output relation between 4 control inputs to a preselected set of a few sensor (say, 5 or 
7) locations These sensor locations have been chosen from an optimal set of 23 sensor 
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locations obtained in [21] following a procedure involving Fisher Information Matrix and 
Effective Independence Distribution vector The optimal 23 sensor locations are as shown 
in fig (3 3) 


1 3 Organisation Of The Thesis 

The chapters m the thesis are organized as follows Chapter 2 presents a description of 
equations of motion of the helicopters, followed by an introduction to the control law 
design The algorithm for flexible/rigid body mode decomposition and the algorithm for 
model order reduction are described This is followed by design of the controller, by 
factorisation theory using internal model principle for disturbance rejection 
Chapter 3 presents the results and discussion of the present study Chapter 4 summarizes 
the results and gives concluding remarks on the present study 
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Chapter 2 

Design And Development of Control 
Laws For Vibration Minimisation 

2 1 Introduction 

For the design of controller for vibration reduction of helicopter at the disturbance fre 
quency the necessary equations of motion describing the system are obtained from [21] 
As the initially given system is of high order [52 order state space system] a suitable 
reduced order model approximating the original system is obtained This necessitates 
the decomposition of the original system into two subsystems The flexible mode subsys 
tern and the rigid body mode subsystem The flexible mode subsystem is then reduced, 
followed by formulating and implementing a strategy for disturbance rejection employing 
the well known internal model principle The control force used for vibration reduction 
is then calculated for the disturbance frequency 

2 2 Description Of The Coupled Gearbox- Fuselage 
System 

For a complete description of the coupled gearbox fuselage system for vibration analysis, 
we need equations of motion describing rigid body and flexible modes of the helicopter 
The equations of motion required are as follws The details are provided in [21] 

1 The rigid body equations of motion of the gearbox It includes 
• Translational equations of motion 
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• Rotational equations of motion comprising of pitch and roll equations 

2 The rigid body equations of motion of fuselage They are as above i e 

• Translational equations of motion 

• Rotational equations of motion comprising of pitch and roll equations 

3 The equations of motion of the elastic modes of the fuselage 

These equations constitute the complete set of linear coupled differential equations used 
for vibration analysis 

In state space form these equation can be written as 


{«} = Ml?} + Ml?} + {/} 


(2 1 ) 


where 


{y} = [C]{g} 


{91} = 




Zgb 

OyGB 

9xGB 

Zp 

OyP 

OxF 

Vl 


Qi 

Q2 


VNm 


{ Q 2 } = {91} 


(6+Wm)xl 


(2 2) 

(2 3 ) 


(2 4 ) 


6 



0(64-Nm)il 


{/} 


EM 

mGB 

0 


0 


(2 5) 


^ 2{6+Nm)xl 




F} 


F! 


(2 6 ) 


NCxl 


where Nm is the number of flexible modes of the fuselage and NC is the number of control 
inputs 

Details of matrices [A] [B] and [C] are provided in the Appendix of [21] [A] being a 
52x52 system matrix [B] a 52x4 control matrix and [C] a 65x52 output matrix 


2 3 Approach to Control Law Design 


Introduction 

The given 52 order system comprises of 46 flexible and 6 rigid body modes or zero 
frequency modes The eigenvalues of the full system are as given in Table(3 9) The 
eigenvalues close to, or almost equal to zero, constitute the rigid body modes 
The given system is decomposed into parallel systems comprising of flexible and rigid 
body mode subsystems The state space model of the decomposed system is given as 



(2 7) 


{?/} = [^7i 



y ^2 


\ 


(2 8 ) 


where (Ai,Si,Ci) is the 46 order flexible subsystem and (A 2 ,B 2 )C 2 ) is the 6th order 
rigid body system The eigenvalues of A\ and A 2 are given in table(3 10) and table(3 11) 
The eigen values of Ai are all in the left half complex plane and hence it is a stable system 
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and eigenvalues of A 2 are all nearly zero with one of the values, shown as negative value 
which IS due to numerical computational error 

For vibration control only the flexible modes are considered which have natural frequency 
close to disturbance frequency The rigid body system is almost uncontrollablle and hence 
IS not considered for controller design 


2 3 1 Algorithm for fiexible/rigid body mode system decompo- 
sition 

This decomposition requires flnding a state co ordinate transformation which results in 
matrix A being block diagonal each of the diagonal blocks having the required eigenval 
ucs See[15,22] 

This problem is solved as 


• first reducing the matrix A to schur form 

• secondly reordering the unordered schur form so that the eigenvalues (diagonal 
elements) are grouped as required 

• thirdly reducing the remaining off diagonal blocks to zero without disturbing the 
diagonal blocks 


The schur form of A matrix is obtained by finding a unitary similarity transformation 
say S, using the QR algorithm Refer [22] 

The matrix A is reduced to its schur form as 


A, = S^AS = 


All Ai2 
0 A22 


where Auand A 22 are upper triangular matrices and their dimensions are determined by 
the required state dimensions of the parallel subsystems After this transformation the 
eigenvalues are distributed randomly between Auand A 22 A transformation matrix say 
G, which distributes them as needed is obtained by means of simple bubble sorting via 
givens rotation, which interchanges a pair of adjacent diagonal elements of S^AS [15] 
After this, we have 
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T = G^AsG = 


^11 21,2 
0 r22 

where Tn and T 22 have the 46 flexible mode and 6 rigid body mode eigen values respec 
tively 

A further state coordinate transformation to reduce T 12 to zero is found This is done 
using a transformation matrix of the form, 

L ° ^ 

in which dimensions of the blocks are the same as T 
As, 


Then 


Y-^TY = T = 


Tn 0 
0 T22 


TY = YT 


Which leads to the Sylvester equation, 


Tn^u — Y12T22 = —Ti 2 


which IS solved foryi 2 Refer [22] 

The final coordinate trasfermation required to decompose the original system [A,B,C D] is 


R = GY 


so that the decomposition is obtained by simply partitioning the matrices of the realise 
tion [R-^AR,R-^B,CR,D] 

2 3 2 Algorithm for Model Order Reduction 

The stable 46th order flexible mode subsystem is now considered for order reduction 
The mam objective of order reduction is to. 
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• Find a reduced order model that reasonably approximates the frequency response 
of the original system upto 4/rev 

• Secondly the controller design for reduced system should effectively stabilise and 
control the full order system as well 


Balanced realisation being one of the best methods for reducing the order of flexible 
space structures, having closely placed eigenvalues Refer [10] we find the balanced 
coordinates for the given 46 order subsystem 

A balanced realisation [Aj B* Ct, Di] is one for which the controllability and observability 
grammians which are solutions of the equations 

A,Wr + WrAl + BbBj = 0 

AjWi, + WbAh + CfCb = 0 


are equal, that is, 

Wr = Wo = E, a diagonal matrix Refer [8 11,17] 

The state coordinate transformation for this is found using the algorithm in [13] and is 
briefly outlined below 


• Compute cholesky factors of the grammians 

• Let Lr and Lo denote the lower triangular cholesky factors of grammians 
Wr = LrL'^ , Wo = LoL^ 

• Compute singular value decomposition of the product of cholesky factors 
L'^Lr = [/ A 

where U and V are unitary matrices and A, a diagonal matrix 

• Form the balancing transformation 
T = LrVA~^ 
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• For the balanced state space realisation 
A = T-^AT B = T-^B C = CT 
we find that 

T-^WrT-^ = T^WoT = A 

After this, by observing the hankel singular values of the full system as given in Ta 
ble(3 12) the part which contained the lowest singular values i e which are most weakly 
coupled to the inputs and outputs were truncated It was found that an 18th order 
model whose hankel singular values are given m Table(3 13) was sufficient to satisfy the 
objectives as mentioned earlier 

2 3 3 Controller Design 

The reduced order model is used for controller design For this we use the standard 
controller parametrisation of the factorisation theory (Appendix A) We initially consider 
the system with 4 control inputs, One disturbance input and 5 outputs This system is 
represented in standard 4 block form as 
where 


w 


z 

u 

G= 

Gii Gi 2 
G21 G22 


y ^ 




Figure 2 1 The Generalized 4-Block Plant 


w = exogenous input consisting of disturbance signal 
u = control signal 
z = output to be controlled 
y = measured output 

A feedback controller C is attached to the plant to form the feedback system shown 
below 

The standard problem is to find a real, rational proper C to minimize the Hoo norm of 
Tju, the transfer matrix from w to z under the constraint that C stabilise G Refer [24] 
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Figure 2 2 The Standard Block Diagram 


Using factorisation theory as explained m Appendix B, the transfer matrix can be 

written as 

T,, = Ti - T2QT, 

where 

Ti T2 and, T3 

are transfer matrices obtained from the plant and C(s) as given in Appendix B 

For asymptotic disturbance rejection of the sinusoidal signal at 4/rev, we consider the 
internal model principle See [6,15,16,26] According to internal model principle, we can 
obtain asymptotic disturbance rejection by incorporating m the feedback path, a suitable 
reduplicated model of the dynamic structure of disturbance signal The purpose of the 
internal model is to supply closed loop transmission zeroes which cancel the poles of the 
disturbance signal 

To attain this using the model matching criterion, the values of Ti T2 and T3 ma 
trices at s = ±j4 is found such that 
Ti(±j4) - T2(i;j4)g(±j4)T3(ij4) = 0 

The value of Q at s = ±^4 which satisfies above equation is found out They are complex 
conjugates of each other 

The state space matrices for Q are found out by equating each term of Q(d:j4) by 
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transfer functions of the form 


s'^+as+bs 

{s+iy 

where a and b are found by solving the two equations given below 
Solving, for s = j4 

= value of respective term at Q(j4) 
and for s = —ji 

value of respective term at Q( j4) we obtain the values of a and b After the 
transfer matrix of Q is found out the stable state space block observer form of Q was 
determined as given m [12] The Q so obtained was used in obtaining the stabilizing 
controller 


The controller designed for 18th order plant was found to stabilise the 46th order plant 
without introducing any new unstable modes when used along with the full system 

2 3 4 Calculation of Control Forces 

For the calculation of control forces, required to be applied by the actuators at the gear 
box mountings, we consider the state space description of the helicopter system and 
controller as given below 

Xp = Axp + Bup + / 

Plant 

Up — Cxp 4 ” DxLp 


— A^Xc 4“ BqXI^^ 

Controller 

Vc — CqXc “b DqUc 

Up IS the plant output, representing the vibration measurement at 5 sensor location 
for a disturbance frequency f, up is the input to the plant, representing the 4 control 
forces applied Uc and yc are the respective controller inputs and outputs respectively 
Now, Vp = Uc 
and 

Vc --Up 
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Figure 2 3 The Closed Loop System 


Up IS calculated as the open loop vibration(ie, Up — 0) of the system for 4/rev input 
disturbance frequency, f i e, 

Vp = [Cp[^p - i'wI]~^Bp + Dp\ f 
where, ui = A/ rev As Hp = Uc, 

From this, knowing the controller description, the controller output force is calculated 
as 

Pc = [Cc[Ac - iw I]-'^ Be + Del Uc 

Up = Vet gives the control force required for vibration minimisation 

Thus the control laws were formulated using the above given algorithms and the 
stabilising controller designed for vibration minimisation ,using the reduced order model 
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Chapter 3 

RESULTS AND DISCUSSION 


This chapter deals with the presentation of results of control laws that were discussed in 
the previous chapter The results of these studies are presented in the following sections 
In the first section a set of 5 sensor locations , chosen from a optimal set of locations as 
given in Ref [21] is considered The results of the closed loop control scheme and open 
loop control design, ( as m [21]) are analysed The second and third sections, deal with 
arbitrary choices of high and low vibration node locations respectively, for closed loop 
control In the fourth section , a sensor at the tail portion of the fuselage is considered 
along with optimal sensor node locations Finally 7 sensor locations from the optimal 
set IS considered for vibration minimisation study 


3 1 Simulation Results For Various Choices Of Sen- 
sor Locations 

3 11 Optimal Sensor Locations (Node Locations 8, 17, 23, 34 
and 39) 

Closed Loop Control Scheme 

The data for mathematical model of the helicopter system in the state space form was 
taken from Ref [21] For the initial study, optimal sensor locations at nodes 8, 17, 23, 34 
and 39 which have high levels of vibration are chosen The actuators are placed at nodes 
39, 48, 46 and 37 The 
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given 52 order system is decomposed into 46 and 6th order subsystems as explained 
m sec (2 3 1) The 46th order system considered for controller design is a stable system 
Balanced realisation based model order reduction technique as explained in sec (2 3 2) 
is employed to transform the given system into its balanced state space coordinates The 
hankel singular values of the 46th order system as given in table[3 12] gives an idea of 
the number of modes that are comparitively strongly coupled to the input and output 
By trial, the system was initially truncated to a 10th order model The 10th order 
system sufficiently approximated the frequency response of the original full system But 
after controller design it is seen that the controller could not effectively stabilise the 
full order system This could be because of spillover instability which is caused by the 
inadvertent excitation of the truncated or residual modes and their unwanted contribution 
to the sensed output Refer [1] and [5] The size of the reduced order model is increased 
one at a time till stability was attained An 18th order model is found that satisfactorily 
approximated the frequency response of the full system and also stabilised the full system 
for the chosen 5 sensor locations This can be seen by comparing fig [3 4] and fig[3 5] 
The individual comparison for all the 5 sensors are shown m fig[3 6] to fig[3 10] The 
lower frequency ranges gave a better match as compared to higher frquency ranges as 
the helicopters natural frequencies are all in (1/rev 6/rev) range 

Following the procedure described m chapter 2, the controller designed is of 22nd or 
der The controller designed for reduced order model is incorporated m the feedback loop 
of the full system and the vibratory levels is calculated at all nodes The vibratory levels 
at all the 5 chosen locations after control as compared to the baseline values is shown 
individually from fig [3 11] to fig [3 15] Their dB values are as shown in the table(3 1) 
given below 


The comparison between baseline and controlled vibration level at all the 64 nodes is 
shown in fig[3 16] The node location 0 indicates the sensor placed at the centre of gravity 
of the gearbox The peak vibration g level in the baseline configuration of the helicopter 
IS 0 2696g at node location 33 and lowest value of vibration is at node 57 of 0 00017g 
The closed loop design, gives maximum vibration at node location 33 of 0 0364g which is 
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Node No 

Baseline 

dB value 

Closed loop 

dB value 

8 

24 5398 

37 3724 

17 

14 1060 

128 157 

23 

21 1613 

31 1831 

34 

13 3981 

29 6488 

39 

21 7099 

33 4949 


Table 3 1 Comparison of Baseline and Controlled Vibratory Levels 


lower than peak baseline value Vibration in the tail portion of the fuselage ( i e nodes 
above 50) is seen to be more after closed loop design For example at node location 64, 
there is a increase in the vibration g level from the baseline value of 0 0032g to 0 0145g 
This could be because there is no sensor chosen from the tail region Vibration at the 
gearbox i e , node 0 is seen to reduce from the baseline value of 0 0598g to 0 01864g 
The control forces are calculated as per the procedure given in sec(2 3 4) The magni 
tude and phase angles of forces applied at each of the actuator nodes is given in table(3 2) 


Node No 

Magnitude 

(10* 10^) 

Phase angle 

(deg) 

39 

19218 

154 19 

48 

2 9892 

-154 53 

46 

0 8335 

52 50 

37 

0 5604 

45 77 


Table 3 2 Control Forces for Optimal Sensor Locations (Closed Loop) 


Open Loop Control Scheme 

An open loop control formulation is also done for the same set of five sensor locations 


17 
















from the optimal set, as per the design given in [21] After open loop control the vi 
bration is seen to reduce significantly at all nodes, other than the gearbox where there 
IS a increase from baseline value of 0 0598g to 0 0624g The comparison between open 
loop- controlled values and the baseline values is shown in fig(3 17) The maximum level 
of vibration after control is at node location 42 of 0 00054g which is quite lower than the 
peak baseline value of 0 2696g at node 33 

The control forces required are much lower than that for closed loop control scheme The 
magnitude and phase angles of the control forces are given below in table(3 3) 


Node No 

Magnitude 

Phase angle 

(deg) 

39 

3 7521 

179 21 

48 

3 7567 

179 21 

46 

3 7558 

179 48 

37 

3 7580 

179 35 


Table 3 3 Control Forces for Optimal Sensor Location (Open Loop) 


3 12 Arbitrary Choice Of High Vibration Nodes (Node Loca- 
tions 12, 16, 32, 33, and 38) 

For the next study another set of 5 arbitrary sensor locations at 12, 16, 32, 33, 38 are 
chosen They are high vibratory level node locations but not optimal as per Ref [21] 
A similar procedure as explained m earlier section is carried out A 20th order model 
was found to sufficiently approximate the full system This gives a 24th order controller 
The comparison between baseline and controlled vibration level for the the full system is 
shown m fig(3 18) The maximum vibration is observed at node location 17 of 0 0326g 
which is lesser than the peak baseline level of 0 2696g at node 33 In the tail region ,there 
IS seen to be an increase in the level of vibration For example, at node location 64 ,the 
vibration increases from baseline value of 0 0032g to 0 017808g In the gearbox, there is 
reduction in vibration from a baseline level of 0 0598g to 0 01359g after control Thus the 
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high vibratory nodes are seen to be effective for vibration reduction The control forces 
are almost similar to forces for optimal sensor locations The magnitude and phase angle 
of control forces are as given in table(3 4) 


Node No 

Magnitude 

(10* 10^) 

Phase angle 

(deg) 

39 

1 1754 

150 54 

48 

1 4224 

137 08 

46 

1 4095 

70 38 

37 

0 9981 

63 31 


Table 3 4 Control Forces for Arbitrary High Vibration Nodes (Closed Loop) 


3 13 Arbitrary Choice Of Low Vibration Nodes (Node loca- 
tions 6, 20, 28, 44, and 64) 

A similar study was carried out for another set of arbitrary 5 sensor locations having low 
vibratory levels Here a 30th order model was needed to sufficiently approximate the full 
order system The controller is of 34th order 

The comparison betwen baseline and controlled vibration level for the full system is 
shown m fig(3 19) The maximum vibration is observed at node 33 of 0 2205g which is 
close to the peak baseline value of 0 2696g At the majority of the nodes the vibration 
after control is seen to be more than the baseline values This may be because the 
sensor measurements of the low vibratory nodes might have been inaccurate resulting in 
inaccurate calculation of control forces There is an increase in the gearbox vibratory 
level from the baseline value of 0 0598g to 0 081 8g 

It IS not safe to choose low vibratory level nodes as sensor location points for vibra 
tion reduction The magnitude and phase of the control forces are as given in table(3 5) 


19 







Node No 

Magnitude 

Phase angle 

(deg) 

39 

43 052 

3 86 

48 

36 73 

-1 96 

46 

30 47 

7 30 

37 

85 52 

010 


Table 3 5 Control Forces for Arbitrary Low Vibration Nodes (Closed Loop) 

3 14 Optimal Sensor Locations With A Tail Sensor (Node Lo- 
cations 17, 23, 34, 39, and 64) 

For this study, the same initial set of optimal sensor locations were chosen with one of 
the nodes 8, replaced by a sensor chosen at the tail of the fuselage i e node 64 A 22nd 
order model approximates sufficiently the full order system 

The comparison between baseline and controlled vibration level for the full system 
is shown in fig (3 20) The maximum vibration is observed at node 38 of 0 02743g 
Vibration is seen to be lowered at all the nodes Specially in the tail region , there is 
seen to be much reduction as compared to earlier choices For example at node location 
55, there is seen to be reduction from the baseline value of 0 00367g to 0 00041g This 
could be because of the presence of a sensor at the tail region along with other high 
vibration nodes 

The control forces calculated here are lower than the 1st and 2nd choices of sensor 
locations They are as shown in table(3 6) 

3 15 Seven Optimal Sensor Locations (Node Locations 8, 17, 
23, 34, 39, 17, and 28) 

Closed Loop Control Scheme 

For this study a set of 7 sensor locations from the optimal set were chosen and the 
controller designed as per the earlier procedure A 24th order system approximates the 
full system 
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Node No 

Magnitude 

(1 O + IO^) 

Phase angle 

(deg) 

39 

0 8474 

165 85 

48 

19157 

165 17 

46 

01161 

11518 

37 

01607 

175 32 


Table 3 6 Control Forces for Optimal Sensor and Tail Sensor Locations (Closed Loop) 


The comparison between baseline and controlled vibration levels is shown in fig(3 21) 
The peak baseline vibration level is 0 2697g at node location 33 The peak controlled 
response is 0 0379g at node location 62 At the gearbox, there is a reduction in vibration 
from the baseline value of 0 05988g to 0 0218g In the tail region for example at node 
64 there is an increase in vibration from the baseline value of 0 0032 Ig to 0 03749g 
The control forces are as shown in table(3 7) The magnitude of control force is seen 
to have been lowered in the case of 7 sensors by a factor of 10 as compared to optimal 5 
sensor locations 


Node No 

Magnitude 
(1 0 * 102) 

Phase angle 

(deg) 

39 

3 4388 

11 57 

48 

5 3327 

12 83 

46 

0 5525 

84 96 

37 

0 4289 

59 77 


Table 3 7 Control Forces for 7 Optimal Sensor Locations (Closed Loop) 


Open Loop Control Scheme 

An open loop scheme formulation is again done for the same set of 7 optimal sensor 
locations, as per the design given in [21] There is a significant reduction m vibration 
at all nodes except at the gearbox after control The comparison between baseline and 
open loop controlled vibratory levels is shown in fig (3 22) At the gearbox i e ,node 0, 
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there is an increase in vibration from baseline value of 0 0598g to 0 06249g Other than 
gearbox the peak vibration is at node location 42 of 0 00028g which is quite lower than 
the baseline value The control forces required are similar to that required for the open 
loop formulation for 5 optimal sensor location The magnitude and phase angles of the 
control forces are as given in table (3 8) 


Node No 

Magnitude 

Phase angle 

(deg) 

39 

3 7553 

-179 47 

48 

3 7553 

-179 69 

46 

3 7561 

179 47 

37 

3 7567 

179 44 


Table 3 8 Control Forces for 7 Optimal Sensor Locations (Open Loop) 

Thus for various choices of the sensor locations, the control laws designed as per the 
algorithms explained in Chapter 2 were implemented and the results analysed 
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Table 3 9 Eigenvalues of the Full System 


SNo 

Eigen values 

1 

0 08875102976272 +17 29791108351503i 

2 

0 08875102976272 17 29791108351503i 

3 

0 02795818992943 +11 93083803276192i 

4 

0 02795818992943 11 93083803276192i 

5 

0 02800282165061 +12 30482819641770i 

6 

0 02800282165061 12 30482819641770i 

7 

0 03475416128688 + 7 05896159Q53937i 

8 

0 03475416128688 7 05896159053937i 

9 

0 03069535817620 + 6 13029696397045i 

10 

0 03069535817620 6 13029696397045i 

11 

0 03189076637836 + 5 82440867391609i 

12 

0 03189076637836 5 82440867391609i 

13 

0 03291029354580 + 5 53104714394176i 

14 

0 03291029354580 5 53104714394176i 

15 

0 02645884964234 + 5 44131232662959i 

16 

0 02645884964234 5 44131232662959i 

17 

0 02529607303838 + 4 92342188521718i 

18 

0 02529607303838 4 92342188521718i 

19 

0 02487619896223 + 4 79801229595323i 

20 

0 02487619896223 4 79801229595323i 

21 

0 02290578844185 + 4 39551765863855i 

22 

0 02290578844185 4 39551765863855i 

23 

0 02422156490082 + 4 76665773658356i 

24 

-0 02422156490082 4 76665773658356i 

25 

0 02234539199842 + 4 42907155283225i 

26 

0 02234539199842 4 42907155283225i 


Contd 
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S No Eigen values 

27 0 01984470001224 + 3 474351 73914139i 

28 0 01984470001224 3 47435173914139i 

29 0 02380946872847 + 3 70094095012000i 

30 0 02380946872847 3 70094095012000i 

31 0 00484154758626 + 0 66366775694405i 

32 0 00484154758626 0 66366775694405i 

33 0 00408480616817 + 0 81456693456129i 

34 0 00408480616817 0 81456693456129i 

35 0 01526006344587 + 3 05898620442318i 

36 0 01526006344587 3 05898620442318i 

37 0 01463144712625 + 2 24567321755696i 

38 0 01463144712625 2 24567321755696i 

39 0 01091992240451 + 2 19632306442407i 

40 0 01091992240451 2 19632306442407i 

41 0 01329022673454 + 2 63113636358172i 

42 0 01329022673454 2 631 136363581 72i 

43 0 01238520927406 + 2 38374970141317i 

44 -0 01238520927406 2 38374970141317i 

45 0 01308997918679 -t- 2 60197904001544i 

46 0 01308997918679 2 60197904001544i 

47 0 00000000000000 + 0 00000006580385i 

48 0 00000000000000 0 00000006580385i 

49 0 00000000000000 0 00000003886751i 

50 0 00000000000000 0 000000038867511 

51 0 00000003638915 
0 00000003638915 
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Table 3 10 Eigenvalues of the Flexible Mode System 


SNo 

Eigen values 

1 

0 08875102976282 +17 29791108351510i 

2 

0 08875102976282 17 29791108351510i 

3 

0 02800282165060 +12 30482819641770i 

4 

0 02800282165060 12 30482819641770i 

5 

0 02795818992936 +11 93083803276182i 

6 

0 02795818992936 11 93083803276182i 

7 

0 03475416128689 + 7 05896159053937i 

8 

0 03475416128689 7 05896159053937i 

9 

0 03069535817621 + 6 13029696397046i 

10 

0 03069535817621 6 13029696397046i 

11 

0 03189076637836 + 5 8244086739161 li 

12 

0 03189076637836 5 82440867391611i 

13 

0 03291029354579 + 5 53104714394177i 

14 

0 03291029354579 5 53104714394177i 

15 

0 02645884964233 + 5 44131232662960i 

16 

0 02645884964233 5 44131232662960i 

17 

0 02529607303838 + 4 92342188521718i 

18 

0 02529607303838 4 92342188521718i 

19 

0 02487619896223 + 4 79801229595323i 

20 

0 02487619896223 4 79801229595323i 

21 

0 02422156490082 + 4 76665773658357i 

22 

0 02422156490082 4 76665773658357i 

23 

0 02290578844186 + 4 39551765863862i 

24 

0 02290578844186 4 39551765863862i 


Contd 
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S No Eigenvalues 

25 -0 02234539199843 + 4 42907155283226i 

26 -0 02234539199843 4 429071 55283226i 

27 0 00484154758625 + 0 66366775694405i 

28 0 00484154758625 0 66366775694405i 

29 0 02380946872847 -h 3 70094095012000i 

30 0 02380946872847 3 70094095012000i 

31 0 01984470001225 + 3 47435173914146i 

32 0 01984470001225 3 47435173914146i 

33 0 00408480616817 + 0 81456693456129i 

34 0 00408480616817 0 814566934561 29i 

35 0 01526006344587 + 3 05898620442318i 

36 0 01526006344587 3 05898620442318i 

37 0 01463144712629 + 2 24567321755691i 

38 -0 01463144712629 2 24567321755691i 

39 0 01091992240451 + 2 19632306442407i 

40 0 01091992240451 2 19632306442407i 

41 0 01329022673455 -t- 2 63113636358170i 

42 0 01329022673455 2 63113636358170i 

43 0 01238520927406 + 2 38374970141317i 

44 0 01238520927406 2 38374970141317i 

45 0 01308997918679 4- 2 60197904001544i 

46 0 01308997918679 2 60 197904001 544i 
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Table 3 11 Eigenvalues of the Rigid Body Mode System 


SNo 

Eigenvalues 

1 

0 00000000000004 + 0 000000072345371 

2 

0 00000000000004 0 000000072345371 

3 

0 00000000000000 + 0 000000016367331 

4 

0 00000000000000 0 000000016367331 

5 

0 00000002792489 

6 

-0 00000002792489 
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Table 3 12 Hankel Singular Values of the Flexible Mode System 


SNo 

Hankel Singular values 

1 

0 50020032528298 

2 

0 49787295252775 

3 

0 39792717752109 

4 

0 39404491912856 

5 

0 35759304563833 

6 

0 35421141956328 

7 

0 34479476287702 

8 

0 34130683900409 

9 

0 29945680703380 

10 

0 29678280169811 

11 

0 29173934551814 

12 

0 29043391374987 

13 

0 25459758631574 

14 

0 25194694270086 

15 

0 22468950975344 

16 

0 22362607838112 

17 

0 17680269962824 

18 

0 17486793266184 

19 

0 12996833234765 

20 

0 12909073516741 


Contd 


28 








29 





Table 3 13 Hankel Singular Values of the Reduced Order System 


SNo 

Hankel Singular values 

1 

0 50008795899589 

2 

0 49775726294606 

3 

0 39792699404819 

4 

0 39404315404613 

5 

0 35768824775163 

6 

0 35419377218320 

7 

0 34477363921226 

8 

0 34128237140741 

9 

0 29205071328056 

10 

0 28869493786649 

11 

0 25465993779617 

12 

0 25208554185246 

13 

0 17677949506505 

14 

0 17482638374724 

15 

0 12314452287337 

16 

0 12191603548932 

17 

0 09301425910005 

18 

0 09167333022519 
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Figure 3 1 Coupled Gearbox Fuselage Dynamic Model 
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Figure 3 2 Baseline Finite Element Model of Helicopter Fuselage 
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Reduced Order System 
Pull Order System 

Figure 3 7 Comparison of Frequency Response of the Full and Reduced Order Model 
(Node No 17) 
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FREQUENCY (Aev) 


Reduced Order System 
Full Order System 

Figure 3 8 Comparison of Frequency Response of the Full and Reduced Order Model 
(Node No 23) 
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Figure 3 9 Comparison of 
(Node No 34) 
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- Reduced Order System 
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squency Response of the Full and Reduced Order Model 
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Reduced Order System 
Full Order System 

Figure 3 10 Comparison of Frequency Response of the Full and Reduced Order Model 
(Node No 39) 
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- Uncontrolled System 
Controlled System 

Figure 3 11 Frequency Response of the Uncontrolled and Closed Loop Controlled System 
(Node No 8) 
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Uncontrolled System 
Controlled System 

Figure 3 12 Frequency Response of the Uncontrolled and Closed Loop Controlled System 
(Node No 17) 
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Uncontrolled System 
Controlled System 

Figure 3 13 Frequency Response of the Uncontrolled and Closed Loop Controlled 
(Node No 23) 
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Figure 3 14 Frequency Response of the Uncontrolled and Closed Loop Controlled System 
(Node No 34) 
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Figiue 3 16 Comparison of Baseline and Closed Loop Controlled Vibratory Levels (Node 
Locations 8, 17, 23, 34, 39) 
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Figure 3 17 Comparison of Baseline and Open Loop Controlled Vibratory Levels (Node 
Locations 8, 17, 23, 34, 39) 
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Figure 3 18 Comparison of Baseline and Closed Loop Controlled Vibratory Levels (Node 
Locations 12, 16, 32, 33, 38) 
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Figure 3 19 Comparison of Baseline and Closed Loop Controlled Vibratory Levels (Node 
Locations 6, 20, 28, 44, 64) 
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Figure 3 20 Comparison of Baseline and Closed Loop Controlled Vibratory Levels (Node 
Locations 17, 23, 34, 39, 64) 
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Figure 3 21 Comparison of Baseline and Closed Loop Controlled Vibratory Levels (Node 
Locations 8, 17, 23, 34, 39, 15, 28) 
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Chapter 4 

Concluding Remarks 


In this thesis, I he formulation and design of a closed loop control scheme for vibration 
mininiis ition of the helicopter has been carried out The full order helicopter system(52 
ordei) wis decomposed to flexible mode subsystem(46 order) and rigid body mode sub 
system(6 order) A reduced order model was found out, that satisfactorily approximated 
the full system, using the technique of balanced realisation based model order reduction 
technique A stabilising controller is designed for the reduced system using factorisation 
theory and internal model principle for disturbance reduction 

For the choice of five optimal sensor locations, the closed loop control gives reduced 
vibration levels as compared to the baseline values, m the fuselage and gearbox In the 
tail region though, there is observed to be an increase m vibration from the baseline 
value This could be because there were no sensors chosen from the tail region On 
having a s( nsor in the tail along with the optimal sensors, the level of vibration m the 
tail portion is seen to be lowered In the open loop control formulation, as compared to 
closed loop design, there is more reduction in vibration, in the fuselage But in the 
gearbox ,it is observed that there is an increase in vibration after open loop control 

An arbitrary choice of high vibration node locations for sensor placement gave good 
vibration reduction as compared to arbitrary low vibration sensor locations Thus, it 
IS observed that sensor locations have a significant influence on the level of vibration 
reduction in fuselage structure 

It IS seen that there is very small variation in the magnitudes and phase angles of the 
control force, in open loop control scheme, on increasing the number of optimal sensor 
locations from five to seven While for closed loop scheme, for seven sensor locations, 
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the control force is seen to be reduced by a factor of 10 as compared with five sensor 
locations The level of vibration reduction is almost similar for five and seven number of 
optimal sensor locations 


54 



APPENDIX A 


Parametrization Of All Stabilising Controllers 


Using f xctorisaton theory Ref [24] we obtain a doubly coprime stable factorization of 
a plant transfer matrix P Using this, the formula for all stabilising compensators which 
can also guaiantee poles of the closed loop system in a specified region of the complex 
plane, if required, is obtained 

Consider the system transfer matrix P represented in a state space model as 


X = Ax + Bu 
y = Cx + Du 


which IS denoted in shorthand as 


P 


^ A B 
D 


= CisI-A)-^B + D 


The doubly coprime factorization of P is defined as follows 
Definition 

A doubly coprime factorization(DCF) of P over a closed region U m the complex plane 
IS a factorization of P as 


P = ND~^ = D'^N 


such that the following properties are satisfied 

1 N, D, N, D have no poles in Q 

2 There ate transfer matrices X, Y, X, Y which have no poles m SI 

3 Y, y aie square and nonsingular 

4 N, D, N, D> X, y, X, Y satisfy the identity 


y X 

-N D 


D -X 
N Y 


= I 
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The procedure for obtaining N, D, iV, D of a given plant P from its state space discription 
IS given below Suppose the given system 


\C D 

IS minimal Select matrices K and H such that Ak - A- BK and Ah = A— HC have 
no eigcn vilues in 0 Then N, D, N, D, X, Y, X Y are given as [24] 


N = 

D = 

N = 

D = 

X = 

Y = 

X = 

Y = 


Ah 

B-l 

C 

D 

Ah 


-C 

^ J 

Ak 

C- 

\ 

DK 

f 

Ak 

B ^ 

^ -K 


Ah 


K ^ 

oj 

Ah 

B- 

[ ^ 

J 

{ Ak 

H '' 

u 

0 ; 


B 



Formula For All Stabilising Controllers 

In the DCF above, if is a compliment of an open connected domain T in the left half 
of the complex plane, then the factors N, D, h?, D, X, Y, X, K are all stable transfer 

functions with all the poles in F 

The well known formula which makes use of the DCF and expresses all the stabilising 
compensators C(3) such that the closed loop system with the transfer matrix from w the 
disturbance input to z, the regulated output has all its poles in F is given as follows Let 
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Q be iny proper rational transfer matrix with no poles in Q, then 


C = { (y - QN)-\X + QD), det{Y -QN)^0} 
= {{X + DQ){Y-NQ)-\ det{Y-NQ)^0} 


Thus we obt un the parametrizaton of all compensators that stabilize a given plant P 
This p xranictnzation is affine in the free parameter Q 
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APPENDIX B 


The Four Block Control Problem 

In the standard block diagram used for four block problem, as shown m fig () is 
(•he exogeneous input, typically consisting of command signals, disturbances and sensor 
noises u is the control signal z is the output to be controlled, its components typically 
being tracking errors, filtered actuator signals etc y is the measured output A feedback 
controller C is attached to the plant to form the feedback system 



For G, the generalised plant given by the four block transfer matrix 


G 


Gii G\2 

Gil G22 


The equations of the feedback system are then given by 


2 = Gnw + G12U 
y = G21IO -h G22'^ 
u = —Gy 

Hence for a well posed system(.e, when det(I+GC)^ 0), the transfer matrnc denoting the 
response z due to exogenous input w is obtained as 
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where 


Tzw = Gn - GuC{I + G22C)~^G3 i 

Standard problem - 

The standard problem is to find a real rational proper C to minimize the norm 
of the transfer matrix from w to z under the consramt that C stabilizes G All the 
requirements like robustness, disturbance rejection, tracking accuracy etc can be captured 
in the transfer matrix between w and z One such procedure using model matching theory 
is outlined below 

Model matching Theory 

The transfer matrix between w and z can be written as [7] 

T.u,=Ti-T2QT3, 


with C(s) given in Appendix A Ti, Ta and T3 are defined as follows Let the standard 
set up of Fig (5 1 ) in Appendix A has state space realization 


G(s) 


/ 




A Bi 
Cl Dn 
C2 D21 


Choose F and H such that Ap = A - BF and Ap 
and T3 aie given by 


Du 
D22 j 

= A~ HC stable Then Tj, r2 


Ti{s) = 


^4 £ \ 

Dn) 


A = 


Ap B2F 
0 Ah 

<■ 


B = 

a - 


Bi 

Bi — HD21 

C'i-A 2 F D12F 
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( 


T2 {s) 


Ap B2 

\Ci-D,2F Di2 


Tz(s) 


Ah Bi - HD21 ^ 

^2 D21 j 
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APPENDIX C 


Model Reduction Via Balanced Realization 


Model reduction is one of the important issues in the dynamic analysis of flexible struc 
tures A model with a large number of degrees of freedom although useful in a static 
analysis, can cause numerical difficulties uncertainties and higher computational costs, 
if used in dynamical analysis of a system 

For state space models, a methodology for deriving reduced order model is provided in 
terms of a system’s realization in balanced coordinates The key computational problem 
there is the calculation of a balancing transformation and the matrices of the balanced 
realization 

A simple way to decrease the order of a model is to delete everything but the controllable 
and observable part The problem is, this part of the model is structurally unstable and 
that any system is in some sense genencally controllable and observable Therefore, 
one would like to be able to measure the ’’degree of controllability and observability” of 
different subspaces of the state space The most controllable and observable part could 
then be used as a low order approximation for the model Ref[17] 

Bruce C Mooie in his paper, Principal Component Analysis m Linear Systems Con 
trollabilitij, Observability and Model Reduction explains this point The controllability 
and observability grammian are used to define measures of controllability and observ- 
ability m certain directions of the state space The grammians are not mvanant under 
coordinate transformations and it is shown that there exists a coordinate systems m 
which the grammians are equal and diagonal The corresponding system representation 
IS called balanced A reduced order model can be obtained from the balanced represen- 
tation by deleting the least controllable and therefore least observable part 



Consider a linear time invariant continuous time system m state space form 

X — Ax + Bx 
y - Cx + Dx 

where A is nxn, B is nxr and C is mxn 

The system is assumed to be controllable and observable, which means that the control 
lability grammians 

Wr= r 

and the observability grammian 

Wo= r 

JO 

arc both nonsingular for any r > 0 The grammians are not invariant under eqmvalance 
transformations on the system These always exists an equivalent system for which the 
grammians are equal and diagonal, le 


Wr Wo = ^ 

Each state variable is, in a sense, equally strongly coupled to the input and the output 
m a balanced realization and this lies at the root of the usefulness of these realizations, 
for model approximations 

Let aj, z = 1, 2, , n be the singular values of the controllability and observability gram 

mians 

Assume that the state variables have been permitted so that a +i < fft, z = 1, 2 , n - 1 

and that the balanced system is partitioned as 




y = (Cl C2) 


+ Du 


where An is kxk 

Let m and 1 x 2 be the L 2 norm functions that drive the state from the ongin to 


2 



(xf (r) and 

(O x’^{t)Y, in the time interval [0, r] 

It was shown by Moore that 

r II ^2 IP dt ^ gfc \\x2{t)\\^ 

Jo II IP dt - II rri(r) |P 

If afc ^ and ui and U 2 have the same norms it follows that 

II ^2(r) ||C|| ii(r) II 

In the other words, the part X 2 of the state is much less affected by the input than the 
part rrj 

Analogously, let yi and 1/2 he the zero input responses from 0^^ 

and 

(0 x'^{t)Y 

Then 

/ II 2/2 (t) 11^ dt< [ II yi{t) IP dt 

Jo Jo 

if CTfc > (Tk^ 1 

and 

ll^i(o) HI ^2(0) II 

This means that thea ;2 part of the state affects the output much less than the Xi part 
Thus w( see that the X 2 part of the state does not affect the input output behaviour 
of the systc m very much if 
CTiii <crk 

The algonl hm to obtain the transformation matrix which gives a balanced realization is 
given m Sec (2 2 2) 

Algorithm for Model Order Reduction 

To obtain an approximate model of the original system a balanced realization of 
the original system is obtained and then simply truncated by discarding those parts 
relating to the state variables which are most weakly coupled to the inputs and outputs 
Ref[13,15] 

A workable algorithm that gives a reduced model of order k after balancing is , 



• Fmd Wc and Wo for the given system (A B, C) 

• Let k = Order of reduced order model 

• Compute cholesky factors of the grammions Wr and Wo 

Wr = ULJ 

and 

Wo^LoLo'^ 

• Compute singular value decomposition of the product of cholesky factors 
LlL, =UaV^ 

where U and V are unitary matrices and A, a diagonalmatrix 

• Find, S = and let 

Ui ~ Matrix containing only the first k columns of U 
Vi — Matrix containing only the first k columns of V 

• Form the balancing transformations, 

T\ = Lo *Ui * S 
T2 = L^*Vi*S 

• Then, the reduced model {Ak,Bk,Ck,D) is, 

Afe = Tf * A * T 2 

Bk = Tl*B 
Ck = C*T2 
Dk = D 


Thus we obtain a balanced, reduced, kth order model 
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